Primordial Bispectrum Information from CMB Polarization 

Daniel BabichQ and Matias ZaldarriagJ3 

Deparment of Physics, Harvard University, MA 02138 and 
Harvard- Smithsonian Center for Astrophysics, 
60 Garden Street, Cambridge, MA 02138 
^ ■ (Dated: February 2, 2008. To be submitted to Phys. Rev. D.) 

8 ■ Abstract 

After the precise observations of the Cosmic Microwave Background (CMB) anisotropy power 

^ : 

, spectrum, attention is now being focused on higher order statistics of the CMB anisotropics. 

\ Since linear evolution preserves the statistical properties of the initial conditions, observed non- 

Gaussianity of the CMB will mirror primordial non-Gaussianity. Single field slow-roll inflation 

. robustly predicts negligible non-Gaussianity so an indication of primordial non-Gaussianity will 

^) \ suggest alternative scenarios need to be considered. In this paper we calculate the information 
00 ■ 

' on primordial non-Gaussianities encoded in the polarization of the CMB. After deriving the opti- 
^ ■ 

■ mal weights for a cubic estimator we evaluate the Signal-to-Noise ratio (S/N) of the estimator for 



WMAP, Planck and an ideal cosmic variance limited experiment. We find that when the experiment 



O ' can observe CMB polarization with good sensitivity, the sensitivity to primordial non-Gaussianity 

^' 

^ ■ increases by roughly a factor of two. We also test the weakly non-Gaussian assumption used to 

■ 

derive the optimal weight factor by calculating the degradation factor produced by the gravita- 
• 1—1 . 

^ I tional lensing induced connected four-point function. The physical scales in the radiative transfer 

functions are largely irrelevant for the constraints on the primordial non-Gaussianity. We show 
that the total (S/N)'^ is simply proportional to the number of observed pixels on the sky. 
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I. INTRODUCTION 



Recent advances in observational cosmology have led to unprecedented constraints on 
the parameters of the cosmological model. One of the outstanding goals in the field is 
to determine the mechanism responsible for the seeds that lead to the the structure in 
the observable Universe. The current data is consistent with an initial scale invariant and 
adiabatic spectrum of primordial curvature perturbations which existed outside the horizon 
at recombination. These observations are in agreement with the predictions of single field 
slow-roll inflationary models, but various alternatives scenarios still remain viable. 

In the next decade we will see further advances in observations that will constrain many 
aspects of the primordial seeds. The spectral index of their power spectrum will be measured 
more accurately and over a wider range of scales by a combination of Cosmic Microwave 
Background (CMB) and other Large Scale Structure (LSS) probes. Dedicated CMB po- 
larization instruments will establish whether or not there is a stochastic background of 
gravitational waves as predicted by models where inflation happens at the GUT scale. The 
constraints on Gaussianity, which are the focus of this paper, will also improve significantly. 

The Gaussianity of the primordial fluctuations has received a lot of attention lately mainly 
because in standard single-field slow-roll inflation models the deviations from Gaussianity 
of the perturbations can be fully calculated and are directly related to the departures from 
scale invariance and thus predicted to be very small Q], most probably unobservable in the 
CMB. Thus constraints on the Gaussianity can help distinguish simple inflation models from 
the various alternatives. Most of the alternatives to single-field slow-roll inflation solve the 
standard cosmological problems in the same way, by invoking a period of accelerated expan- 

n 

sion in the very early universe (see [2] for a counter example). They differ however in the 
characteristics of the produced perturbations. Just as in slow-roll models the perturbations 
arise from quantum fluctuations during inflation but the detailed physics is not the same. 
In some of the models the dynamics of the inflaton field is fundamentally changed by the 
presence of higher derivative terms in the Lagrangian which can even lead to an inflaton field 
that is not rolling slowly 0, S, 0] . Another possibility are models in which fluctuations 
in another field different from the inflaton are responsible for the adiabatic fluctuations we 
observe today 3, Ol • 

These alternative models usually make distinctive predictions about the shape of the 
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spectrum of primordial perturbations, the amplitude of the gravity wave background and 
the departures from Gaussianity as described, for example, by the 3-point function. For 
scale invariant perturbations the three point function in Fourier space or bispectrum, is 
effectively a full function of two variables, thus it could contain a wealth of information 
about the primordial seeds. The different alternatives to slow-roll inflation not only predict 
different levels of non-Gaussianities but also different shapes for the bispectrum as a function 
of triangle configuration. Moreover one can make a very definite and model independent 
statement about the shape of the three-point function in the so called collapsed limit. The 
collapsed limit corresponds to a three-point function where one of the Fourier modes has a 
much longer wavelength than the other two. In that limit the three-point function should 
;o to zero, unless more than one degree of freedom is dynamically relevant during inflation 
l]. Other consistency relations can be obtained involving the three-point function and the 
amplitude and slope of tensor perturbations [lol |. 

Thus a detailed measurement of the three-point function could provide very interesting 
information on the mechanism responsible for generating the primordial curvature pertur- 
bations. In this work we will use the so called "Local Model" for the non-Gaussianities and 
specify the amplitude of non-Gaussianity by the parameter /^tl. In this model the gravi- 
tational potential, ^{x), can be expressed in terms of a Gaussian gravitational potential, 
^g{x), as 

^X) = %ix) + fNLl^lix) - i^lix))]. (1) 

The bispectrum in this model can be written as 

($(fci)$(fc2)$(fc3)) = 2fNL{27r)''6^'\ki + fc2 + ks)[P{ki)P{k2) + eye], (2) 

where P{k) is the power spectrum. 

In general the primordial bispectrum cannot be written as in Eq. (PJ and so constraints on 
ihe /tvl parameter do not automatically apply to all models of primordial non-Gaussianity 



Nevertheless, we can list approximate values of Jnl that are expected in the various 
alternatives to single field slow-roll inflation in order to estimate the detectability of these 
various models. For comparison in slow-roll inflation models one expects an /tvl ~ 0.05 Q|. 
In models where higher derivative operators are important for the dynamics of the field one 
can expect results ranging from /tvl ~ 0.1 in models where high derivative operators are 
suppressed by a low UV cut-off ^i] to /at^ ~ 100 in models based on the Dirac-Born-Infeld 



effective action. A recent idea called ghost inffation, where inflation occurs in a background 
that has a constant rate of change instead of a constant background value, can also give 
/■/yr ~ 100 yi. The effect of additional light fields on the efficiency of reheating can lead 

/a^l > 5 12]. Another example is the curvaton model where isocurvature fluctuations in a 
second light scalar field during inflation generate adiabatic fluctuations after the inflationary 
epoch is completed, these can cause large non-Gaussianity, Jnl > 10 ||8|. 

There are two additional sources of non-Gaussianity that although not primordial in origin 
could be observed first by future experiments. The first and perhaps the most important for 
observations is secondary anisotropics such as gravitational lensing, the thermal and kinetic 
Sunyaev-Zeldovich effects and the effects of a patchy reionization. The second source are 
non-Gaussianities related to the non-linear nature of General Relativity. The expectation 



13 



14|). although there is no full 



is that this latter source will lead to f^i ~ 1 (eg- 
calculation of the temperature and polarization anisotropics of the CMB to second order 
in the the curvature perturbations. In results were presented for the collapsed limit, 
where the calculation simplifies enormously. It that case f^L ~ 0.7 was obtained although 
the shape dependance of the bispectrum was different from that of the /tvl model. 

While some of the afforementioned models of the early universe are speculative, the second 
order single field, slow-roll inflation and non-linear radiative transfer calculations should be 
taken as definite predictions of standard cosmology. It is interesting to determine how close 
we are to observing any of these effects. Currently the best constraints, which come from 
the Wilkinson Microwave Anisotropy Probe (WMAP) Q, imply -58 < /nl < 134 (95% 
C.L.) [17|. 

The theoretical ability of CMB temperature maps to constrain f^i has been determined 
for a COBE normalized Einstein-de Sitter (EdS) model It was shown that the minimum 
Jnl detectable by WMAP is 20, Planck is 5 and an ideal experiment (no noise and infinites- 
imal beam width) is 3. The limiting factor in the case of an ideal experiment was taken to 
be the effect of gravitational lensing, which increases the estimator noise without affecting 
the CMB bispectrum signal. Gravitational lensing adds to the cosmic variance portion of 
the noise because the particular realization of the large scale structure that lenses the CMB 
is also a priori unknown. Other secondary sources of anisotropics will have a similar effect. 

Since it appears that we are close to being able to detect some of the interesting non- 
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Gaussian effects mentioned above, it is important to explore other sources of information 
beyond tfie CMB temperature fluctuations to see if the minimum detectable value is lowered 
when the new information is included. The CMB is linearly polarized and the E type 
polarization is sensitive to the primordial curvature fluctuations. Therefore in this paper 
we will ask how much stronger would the constraints on f^^ be if information from the E 
polarization is included. Since B polarization cannot be directly generated by the scalar 
primordial curvature fluctuations we will ignore it in this paper. 

As we mentioned above, there are two ways in which secondary anisotropics can com- 
plicate and degrade our ability to constrain primordial non-Gaussianity. If the source of 
secondary anisotropics being considered produces a three-point function it might bias the 
estimator of primordial non-Gaussianity. In practice, a model of primordial non-Gaussianity 
is assumed and a reduced bispectrum template is calculated. The distinct shape of the pri- 
mordial reduced bispectrum may allow it to be distinguished from the other bispectra |l8[ , 
assuming it can be detected with a sufficiently high signal to noise ratio. In addition sec- 
ondaries add fluctuations which increase the variance estimators of the three-point function 
without contributing to the signal. In fact the variance of the estimator is related to the 
six-point function of the temperature field. The six-point function can be expressed in terms 
of its unconnected Gaussian contribution: permutations of three two-point functions, the 
product of a two and a connected four-point functions and a connected six-point function. 
Previous work on this subject has ignored any non-Gaussian contribution to the six-point 
function. 

The secondary anisotropy that will be considered in detail is gravitational lensing. Grav- 
itational lensing does not produces a three-point function, so it will not bias our estimator 
of the primordial bispectrum. However it will produce corrections to the two-point function 
and create four and six-point functions even if the CMB anisotropics are perfectly Gaussian 
[ifll l2n| . Also, lensing will create small scale power in the CMB two-point function which 
acts like noise in the analysis of primordial non-Gaussianity. All information about Jnl is 
eliminated on scales where this effect dominates; this limits the minimum fMr that an ideal 

n 

experiment can detect to be f^i ~ 3 for the EdS model |18]. Fortunately the graviational 
lensing effect in a EdS cosmolog y is much larger than the gravitational lensing effect in a 
AC DM concordance cosmology j2l|, so gravitational lensing will be less important in our 
calculations. For reference the lensed and unlensed CMB power spectra are plotted in Fig. 
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FIG. 1: The lensed (dashed red) and unlensed (sohd black) Ci for the concordance AC DM 
cosmology; the temperature power spectra are on the left and the E polarization power spectra 
are on the right. 

^ Lensing increases the power in the temperature (polarization) anisotropies by a factor 
of two at / ~ 4200 (/ ~ 5000). However we will find that the four-point function starts 
reducing the signal to noise of the bispectrum on significantly larger scales. 

In this work we will ignore other secondary anisotropies such as the thermal and ki- 
netic Sunyaev-Zeldovich effects and the Ostriker-Vishniac effect. The dominant secondary 
anisotropy, the thermal Sunyaev-Zeldovich effect, has a characteristic spectral shape that 
will allow it to be separated. In general the physics of these secondary anisotropies requires a 
non-linear analysis using hydrodynamical simulations, we will leave a detailed analysis of the 
effects of these secondary anisotropies to a future work [2^. These secondary anisotropies 
also produce polarization in the CMB, fortunately the amplitude of the polarization sec- 
ondary anisotropies is much lower than the corresponding temperature anisotropies and will 
be less of a problem. 

Finally, although the E and B polarization modes most directly correspond to the pri- 
mordial curvature fluctuations and gravity waves, they are not directly measured in CMB 
experiments. The Stokes' parameters Q and U are measured and then decomposed into E 
and B. While this decomposition is perfectly well defined for a noise free experiment ob- 
serving the full sky, ambiguities arises in practical experimental situations [2^. Fortunately 
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the majority of the information on /jvi comes from small scales, where this ambiguity is 
less of a problem if we assume the beam is oversampled in order to reduce the effects of 



power aliasing that can also mix E and B modes [2J, |25[ . In this paper we will ignore these 
complications in the decomposition of the Stokes' parameters. 

The paper is organized as follows: in section II we derive the optimal estimator when 
polarization information is included in the non-Gaussianity measurements, provide numeri- 
cal results for the improved constraints on /tvl and quantify the reduction in the S/N due 
to the connected four-point function contribution to the noise. In section III we analyti- 
cally reproduce our calculations for toy models that do not include the effect of radiative 
transfer or the curvature of the sky, but develop intuition about our numerical results. We 
will conclude in section IV with a discussion of our results. In this paper we assume the 
standard ACDM cosmology with = 0.73, fife = 0.044, fi^ = 0.226, ag = 0.8 and Hq = 72 
km Mpc~^ with a scale invariant primordial power spectrum and normalized to the 1-yr 



WMAP data 







27l |. Also we do not include the effects of reionization nor the late-time 



integrated Sachs- Wolfe effect. The late-time integrated Sachs- Wolfe effect only affects large 
scales which do not contribute much to the total signal. 



II. POLARIZATION 



A. Optimal Estimator 

The Signal-to-Noise ratio (S/N) defined in can be generalized to include polarization 
information by finding the optimal weight functions for a cubic estimator. First we form the 
estimator of the CMB temperature and polarization bispectrum signal as 

i,j,k {l,m} 

where the indicies k run over T and E, {/, m} indicates all three /j, rrii and W^f^y is the 
weight function we will optimize. The expectation value of the estimator is 

i,j,k {l,m.} 
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where following the standard notation in the literature we separate {ai^rniai^rn20'hmz) into 
the reduced bispectrum h^i^i.^ and the Gaunt Integral, 



^-7^-3 ^ ; (2/l + 1) (2/2 + 1) (2/3 +1) Ml ^2 /3 h h h , 

* ^'^ loOO/lmi m2 ms 

which characterizes the angular dependece of the bispectrum. The sum over z, j, k includes all 
eight possible bispectra {TTT,TTE,TET, ETT,TEE, ETE, EET, EEE}. In the weakly 
non-Gaussian limit we neglect the contribution of the primordial bispectrum to the variance 
of S, which becomes 

= E E E E [cov]j^5;-, (6) 

i,j,k p,q,r {l,m} {t,s} 

here we have defined a covariance matrix between the eight possible bispectra for each value 
of {l,m} and {t,s} as 

After restricting the indices such that h < h < h and ti < t2 < t^, Eq. ((Zj) is evalulated 
using Wick's Theorem in terms of Cf, Cf' or , which are the temperature, E-mode 
polarization and cross correlation power spectra respectively. It is necessary to include 
permutation factors when some of the Z's are equal; we include a factor of two if two Vs are 
equal and six when all three Vs are equal We include instrument noise in the standard 



fashion 29| and adopt parameters values (beam width and pixel noise) that are relevant for 
WMAP and Planck 

Once we have chosen the ordering of the multiples indices, the evaluation of Eq. (|7j) 
produces Kronecker 5's which allow us to rewrite the variance as 

(^') = E E E ^iSi^S i^^Sr- (8) 

i,j,k p,q,T {l,m} 

Defining the (S/N)"^ as (5*)^/ (A^^) we find the optimal weights by maximimizing this ratio: 

6{S) _ (S) 6{N^) 
6W {N^) 6W ' 

where 



(9) 



6{S) 



and 



6W 



2E^s[cHs:ir- (11) 



p,q,r 

It is clear that Eq. © is satisfied when we choose 



= E ^^i^'^'^i. \cov-x:Sr- (12) 

ijk 

Now defining the quadratic form 

{Q^C-^Qhi.i. = EE^^;^i[Cov-^]S:|r?^^L> (13) 

p,q,r 

where q is a vector that contains all the possible bispectra, and using the summation prop- 
erties of the Wigner 3j symbols we find the formula for the S/N when we optimally include 
both temperature and polarization information about the observed CMB, 

S ,2 _ (2/i + l)(2/2 + l)(2/3 + 1) ( h h h 



This is a straightforward generalization of the formula in jisl, where the quadratic form 
in Eq. (HH) is simply replaced by bf^i^iJCi^Ci^Ci^. 



iV^ 2- 4^^ 1 . . . ^ ^)'^'^'3- (14) 

2<Zl <i2 <'3 <'maa: 



B. Reduced Bispectrum 

We follow the notation of in our calculation of polarization and cross correlation 
reduced bispectra. Working to linear order in radiative transfer, the CMB temperature and 
polarization fiuctuations can be expressed as 

«L = 4vrz' j j^^Ai{k)YCMm, (15) 

where j = T,E corresponds respectively to either temperature or polarization, $(k) is the 
primordial curvature fiuctuation and A;(fc) is the radiation transfer function calculated by 
CMBFAST. We include the relevant factor of -y/jj^fji the polarization transfer function. 
The reduced bispectrum, defined in Eq. (jH), can be expressed as 

= / r'dr[(3l{r)(3l{r)al{r) + A\(r)<(r)A^(r) + al{r)(3l{r)(3^^{r)], (16) 

where 

(3i(r) = ^J edkPik)Mkr)Ai{k), (17) 
9 



and 

«K0 = IJ k'dkji{kr)Alik), (18) 

again where i = T or E. In our notation f3i{r),ai{r) is respectively the equivalent of 
6|*"'(r), 6"°'^(r) in the notation of Q|. Defining tq as the present day value of conformal 
time, Tn as the value at decoupling and r^i = tq — tr as the comoving distance to the 
surface of last scattering, the region of integration for r is over the sound horizon (from tq 
to To - 2rji). 

In Fig. 121 we see features familiar from the CMB temperature anisotropy and polarization 
power spectra. The upper level displays the temperature and polarization /5i(r)'s. On large 
scales Pf{r) approaches a constant value independent of the radial distance within the sound 
horizon. The dominant mechanism on these scales is the Sachs- Wolfe effect so the radiation 
transfer function approaches Af{k) = —ji{kr£))/3. Since there is no polarization Sachs- 
Wolfe effect, we see Pf{r) on large scales. As we approach smaller scales we see the 
analogue of the familiar acoustic oscillations. 

C. Numerical Results 

The curve of S/N vs. Imax will give an estimate of the minimum /jvl measurable when 
both temperature and polarization information, as well as their cross-correlations, are in- 
cluded in the analysis. The S/N is linear in Jnl, so we can determine the minimum sta- 
tistically observable f^L by requiring S/N = 1. The CMB is only partially polarized so in 
experiments with relatively large noise it should be easier to measure bispectra with fewer 
E^s than those involving three Es. For example WMAP will measure TTE much better 
then EEE. Therefore it is practical to calculate the change in the constraint on /atl as 
we include bispectra with one additional E. In Fig. IHlthe S/N is forecasted for WMAP, 
Planck excluding the effects of gravitational lensing and in Fig. |3]for an ideal experiment 
(no instrument noise, infinitesimal beam width) also excluding the effects of gravitational 
lensing. 

Using only temperature information WMAP will be able to detect an f^^ of 13.3 and 
Planck an /tvl of 4.7, thus approximately recovering the results of the previous analysis jisl ]. 
Our results are summarized in Table HI When the CMB polarization is measured with good 
sensitivity, it appears that the inclusion of polarization roughly increases the sensitive of 



10 




10* 



10 100 1000 

1 



10* 




FIG. 2: Upper Left: A plot of Pf{r) vs. I for values of r = tq — tr (black, solid line), r = tq — 0.6tji 
(red, dashed line) and r = tq — 1.4r/j (blue, dotted line); Upper Right: same but for Pp{r); Lower 
Left: same, but for aj{r)] Lower Right: same, but for af{r). 

the experiments by a factor of 2. For an Ideal experiment the minimum /^vl is lowered to 
L6, which is close to the predicted size of the corrections due to second order corrections to 
gravitational and hydrodynamical evolution of the CMB. 

In Fig. |3]we see the S/N curve continues to rise with a constant slope for the case of 
an ideal experiment. One might wonder why the physical scales in the radiation transfer 
functions, like the sound horizon or Silk length at the surface of last scattering, do not 
strongly influence the slope of the S/N curve as l^ax increases. The Silk length is defined 
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FIG. 3: All figures are {S/N)fj^^ vs. Imax excluding the effects of gravitational lensing for TTT 
(dot dashed green), EEE (dashed red), TTT+TTE (dashed light blue), TTT,TTE+TEE (dotted 
blue) and all bispectra (solid black). Left: WMAP; Right: Planck. 
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FIG. 4: Plotted is {S/N)fj^^ vs. Imax for an ideal experiment (no instrument noise and infinites- 
imal beam size) without the effects of gravitational lensing for TTT (dot dashed green), EEE 
(dashed red), TTT+TTE (dashed hght blue), TTT,TTE+TEE (dotted blue) and ah bispectra 
(solid black). 
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Exp eriment 


TTT 


EEE 


TTT,TTE 


TTT,TTE,TEE 


All 


WMAP 


13.3 


314 


11.2 


10.9 


10.9 


Planck 


4.7 


8.9 


3.4 


3.0 


2.9 


Ideal 


3.5 


2.6 


2.2 


1.8 


1.6 



TABLE I: Minimum values Jnl detectable with signal to noise ratio of one using stated bispectra 
with WMAP, Planck, and an Ideal experiment. 

as the average distance a photon random walks before recombination. This signals a break 
down of tight coupling between the baryons and the photons that effectively smooths out 
the fluctuations on small scales. One could worry that this effectively introduces an average 
that through the central limit theorem would reduce the level of non-Gaussianity of the 
resulting anisotropies, however the linear growth of S/N with l^ax in Fig- M contradicts 
this intuition. In section III we will use toy models in an attempt to better understand our 
unexpected results. 



D. Gravitational Lensing Corrections 

As we mentioned in the Introduction, the variance of the bispectrum estimator is the six- 
point function which includes contributions from connected two, four and six-point functions 
produced by gravitational lensing. It is straightforward to include the influence of the two- 
point, as the unlensed CMB power spectra are simply replaced with their gravitationally 
lensed counter par ts in the expressions for the variance of the estimator. This approach has 
been taken in The four-point function can become large on small scales |2^, so it 
is important to check how our results change once these contributions are included. In this 
subsection we will investigate the effects of the gravitational lensing connected four-point 
function on the S/N of the estimator we have defined. Since the creation of a connected 
four-point function is the same order in the lensing potential expansion as the creation of 
small scale power in the two-point function we will the treat the two effects together in this 
subsection. 

In principle we could also include the effects of the connected six-point function induced 
by gravitational lensing. The leading order (in the projected lensing potential) contribution 
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to the connected six-point function is fourth order in the lensing potential and thus it should 
be subleading. In addition, to this order in the lensing potential expansion there are extra 
terms in the connected two and four-point functions that should be included. Here we will 
ignore these terms and focus on those second order in the lensing potential. 

We derived the optimal weights assuming the six-point function could be evaluated solely 
in terms of permutations of two-point functions. Once we include the contribution from the 
connected four-point function the variance of our estimator will increase and its S/N will 
decrease. In this subsection we will determine the size of this effect. Similar work has been 
pubhshed on the effects of graviational lensing induced non-Gaussianity on the analysis of 
the B-mode polarization power spectrum 3^. There it was shown that the non-Gaussianity 
reduces the information contained in the B-mode power spectrum. 

For simplicity we will only include CMB temperature information, but the inclusion of 
E polarization is straightforward. Also we will work in the flat-sky approximation which 
is excellent for the relevant small scales. The relationship between the all-sky and fiat-sky 
formalisms is well understood j^^. Here we adopt the conventions 

{a{h)a{l2)) = {27if5^^\h + l2)C{h), (19) 
{a{h)a{l2)a{h)) = {2tx f 5^^\h + h + h)B{hMM)- (20) 

In the fiat-sky formalism our estimator of the three point signal is defined as 

^ = / (^ (£y^ (Sp^(^^' I3)«(li)«(l2)a(l3), (21) 
where W{\i, I2, 13) is the weight function. 

Using the above conventions the expectation value of the estimator is 

= / (2?)2 (fj)2 (21^)2 ^(1^' l3)(2vr)^^(^)(li + h + h)B{\,, 1„ I3), (22) 
and its variance is 

x(a(Ii)a(l2)a(l3)a(i;)a(l^)a(I'3)). (23) 

In the weakly non-Gaussian regime, meaning that the six-point function can be expressed 
in terms of permutations of two-point functions, the estimator variance becomes 

14 



This assumes that there is no strong source of non-Gaussianity. While the primordial non- 
Gaussianity is rather small there is the possibility that secondary anisotropies may cause 
significant non-Gaussianity and must be included in Eq. (j24j) . 

The {S/Ny, defined as {j^Y = {S)'^/{N^), must be maximized by choosing the appro- 
priate weight function, I2, 13). We will find the optimal weight function when in the 
weakly non-Gaussian limit the six-point variance is solely determined by the Gaussian contri- 
butions and characterize the reduction in the S/N of the estimator once the non-Gaussianity 
induced by gravitational lensing is included. Of course we could derive the optimal estimator 
including the connected four-point function in the variance. Indeed if we discover that the 
gravitational lensing four-point function significantly reduces the S/N we should modify our 
weight function. Maximizing the signal-to-noise ratio, in an analogous fashion to above, we 
find that we should choose our weights such that 

WM.U)-i2.n.h,^^ ^^^^^y (25) 
This leads to the standard formula for the S/N in the fiat-sky formalism 

f^/'^)^ - w / . 1. . w^clMilk) ■ '''' 

However there are corrections to the variance in Eq. (j^^ from the connected four-point 
function 

(a(li)a(l2)a(l3)a(l4)) = i2nY5^^\h + I2 + I3 + l4)T(Ii, I2, 13, 14). (27) 
Including these terms, the correction to Eq. becomes 

S{N')= ^ I d\d\d\dXd%d%W{h,h,h)W{\[XX) 

xS^'\h + I'i)C(/i)(5(')(l2 + I3 + 1^ + l^)r(l2, 13, 1^, I3), (28) 

where the factor of nine comes from the cyclic permutation symmetries. The connected 
four-point function induced by gravitational lensing has been previously determined in the 
fiat- sky formalism [igj ]. 

We will determine the effects of gravitational lensing by calculating the reduction in the 
estimator signal to noise when the weights derived by ignoring gravitational lensing are used. 
This means that the weight defined in Eq. (j25j) will always be evaluated using unlensed CMB 
power spectra. However the Gaussian and non-Gaussian six-point functions in Eqs. ()24|) . 
will be evaluated with the lensed CMB power spectra when specified. 
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FIG. 5: The ratio of the {S/N)gl including the various forms of gravitational lensing to the 
{S/N)q excluding the gravitational lensing is shown above. The curves indicate the inclusion in 
the estimator variance of the four-point function, but not the two-point function (dotted red), 
both the four and two-point functions (solid black) and two-point function, but not the four-point 
function (dashed blue). 

In Fig. El we plot the ratio, R = {S/N)gl/ {S/N)o, that determines the effect of gravi- 
tational lensing on our ability to observe primordial non-Gaussianity. We have defined the 
reduced S/N that takes into account gravitational lensing as 



where {N'^)gl includes the gravitational lensing two and four-point functions as specified. 
The solid black curve indicates that the estimator variance includes both the two and four- 
point function created by gravitational lensing. While the dashed blue curve only includes 
the effects of the two-point function and the dotted red curve just the four-point function. 

Figure El shows that gravitational lensing does not affect the ability of WMAP to constrain 
primordial non-Gaussianity, however Planck's theoretical ability will be reduced by nearly 
~ 25%. The figure also shows that including the gravitational lensing correction to the 
two-point function leads only to a minor change in the signal to noise on these scales, the 




(29) 
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leading effect coming from the four point function. In fact if one only includes the lensing 
effect through the two-point function one obtains R > 1 for some Imax which clearly is 
unphysical. Eventually, at high enough / (/ ~ 4000) when the lensed CMB power spectra 
are significantly larger than the unlensed ones, the corrections to the two point function will 
also decrease the signal to noise. 

Recall that we assumed the CMB fluctuations were weakly non-Gaussian when we derived 
the optimal weight functions used to construct our estimator. In a future work we will derive 
weight function including the effects of non-Gaussianity due to gravitational lensing 

III. SCALING FORMULAE FOR THE SIGNAL-TO-NOISE RATIO 

The shape of the S/N vs. Imax curve in Fig. 0] implies that our naive expectations of the 
effects of the photon diffusion are incorrect. We do not see a saturation of the S/N curve 
that we predicted earlier. In order to verify these counterintuitive results we will alter the 
CMBFAST transfer functions by hand to test how sensitive the results are to the actual 
form of the transfer functions. In the first case we explicitly include additional damping by 
multiplying the CMBFAST transfer functions by an exponential Ai{k) —>■ /\i{k)e^'^^/'^°^'^ , 
where the Silk damping scale, k£, = bOO/r^,, is chosen such that the effects appear near 
I ~ 500. In the next case we remove all influence of radiative transfer by choosing, A/(fc) — *■ 
ctjiikro), where a is some constant that will always cancel in the formula for the (S/N)"^, 
Eq. (HH). 

In Fig. ini all the S/N curves are roughly parallel meaning that our radical altering of the 
transfer function only changes the numerical coefficients in the expressions for S/N. The 
functional dependence on l^ax appears to be close to the same for all three examples. This 
suggests that we can understand these scalings by using simple toy models. 

A. No Radiative Transfer 

As a first example we will calculate the S/N in the flat sky approximation and with no 
radiative transfer; therefore we will simply observe the underlying modes restricted to the 
plane of the sky. 
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FIG. 6: {S/N)f^l vs. Imax for TTT of the full calculation including the CMBFAST transfer 
functions (solid black), no transfer function (long dashed blue), the CMBFAST transfer functions 
with additional damping (dashed red). 

Here we adopt the following conventions 

(<l>(fci)$(fc2)) = (27r)35(3)(fci + k2)P{k), (30) 

(<l>(fci)$(fc2)$(fc3)) = (27r)35(3)(fci + k, + k,)B{h, h), (31) 

and maintain the previously specified conventions for the CMB anisotropy Fourier coeffi- 
cients, Eq. Using our assumptions of the flat sky approximation with no radiative 
transfer, the Fourier coefficients of the temperature anisotropies can be expressed as 

^^me''''-S^'\l-kKn), (32) 

where is the distance to the surface of last scattering and fc" is the Fourier wavevector 
parallel to the surface of last scattering. The power spectrum of this model is 

^!^ = 4 = ^^.A^ (33) 

where A is the amplitude of the scale invariant power spectrum P{k) = A/k^. Likewise we 
can find the bispectrum 

B{h,h,h) = '^^^{^ + cyc.). (34) 
18 



Using the standard formula ^'^^^(O) = fsky/i^ and substituting these results in Eq. gives 
(|)2 ^ [suSjA I d%d%d%6^'\h + h + hmiij^, + cycf, (35) 

t/ 12 

and evaluating the above expression we find 

dr = ^Js,yf^,AtaJn^f^. (36) 

n 

The logarithm is typical of scale invariant primordial power spectra [36]. If the primordial 
perturbations were generated by a Poisson process so each point in space was statistically 
independent, the logarithm would be absent and the dependence on Imax would solely be 
t^ax- Equation (jHUjl can be written in a more physical way by relating it to other observables, 

(|)^ = 8(/^,A)^iV,.ln^, (37) 

where Npi^ = fskyl^aax is the number of observed pixels. 



B. Silk Damping: Toy Model 

There remains the question of why physical processes like Silk damping or cancellation 
due to oscillations during the finite width of the last scattering surface do not cause a strong 
change in the slope of S/N curve at high /. First of all it is important to note that there are an 
equal number of transfer functions in the numerator and denominator of Eq. (pUj) . so there is 
a sense that the effects of radiative transfer cancel out. Of course the transfer functions are 
not simple multiplicative factors that can be cancelled. However the saturation of S/N when 
instrument noise or gravitational lensing comes to dominate the Q's, as contrasted with the 
continued growth of S/N in the case of an ideal experiment without gravitational lensing, 
can be understood from this perspective. When gravitational lensing and other secondaries 
or instrument noise dominates the six-point function in Eq. ()26p the convenient cancellation 
cannot occur and we can no longer recover information about primordial non-Gaussianity 
on these scales. In the case of lensing it might be possible to improve the signal to noise 
by using the B type polarization, which on small scales is only generated by lensing, to 
constrain the deflection angle and at least partially "unlens" the observed T and E flelds. 

We will attempt to explore this in the model by including the effects of Silk damping 
by introducing an exponential cutoff to mimic the effects of Silk damping on the radiation 
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transfer function, 

ail) = (27r)2 / -fJl-^k)e"''r^S^^\l - khn)e~''"^'"'"^ , (38) 

where ko is wavevector corresponding to the Silk length below which the radiation transfer 
function is strongly damped. Repeating the above steps we find the power spectrum can be 
formally evaluated in terms of Hypergeometric U-functions as 

C{l) = ^e-^'/^'^U{l/2Aiyil), (39) 

where Id is the 2-d Fourier multiple corresponding to the Silk damping scale as Id = TDko- 
We can make an approximation in order to better understand the effects of Silk damping 
on the CMB power spectrum by cutting off the integral at /c ~ /^d, then 

Cil) = — , (40) 

so when I <^ Id recover Eq. (jHHj) . Likewise we can evaluate the three-point functions 
again in order to facilitate the evaluation of this integral assume that the exponentials cutoff 
the region of integration at fci, ^2 ~ ^d- 

B{IM,) = g%l!e-('?+'^+'i)/^^H / / +cyc.]. (41) 

Then substituting Eq. ^ and Eq. (jH)) into Eq. ^ and assuming I ^ Id 

(|)^ = l^^y^^ J d%d%d%S^'\h + h + ^3)^^^|J^, (42) 
we find that leading term scales as 

(§)' OC fskyfNLAlLx- (43) 

The dependence on l^ax in Eq- is nearly as strong as that in Eq. (jHUj) . This shows 
that we can still expect to recover information about f^i on scales where photon diffusion 
is exponentially damping the transfer functions. In practice, both detector noise, angular 
resolution and secondary anisotropics will limit the smallest scale that could be used. 

C. Physical Arguments 

The strongest feature of Silk damping in our toy model, the exponential damping of the 
CMB power spectrum and the bispectrum, cancel in the expression for the {S/N)'^, Eq. 
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). However the exponential damping is not the only feature caused by this effect, now 
the power spectrum Eq. ()40|) scales like instead of at large /. While the {S/N)'^ 
still scales with Npi^ = fsky^maxi the exact numerical coefficient that determines the slope 
will be reduced. 

We can understand this behavior by considering the contribution of collapsed triangles 
-C ^2 ~ ^3- In this limit the estimator variance, ignoring factors is simply, 

<i,i,-{Ci,+Ni,){Q,+Ni,)\ (44) 

where A*"; represents both the instrument noise and any secondary anisotropy that will de- 
grade our ability to recover the signal of the primary anisotropy. Ignoring numerical factors, 
the reduced bispectrum from Eq. ()16|) can be rewritten as 

fciil ~/7VL j kfdhkldk2kldhAlih)Ai{h)Al{^^^^^ 

(45) 

where we define Cij^izisiki, k2, k^) = f r'^drii^{kir)ii^{k2r)ii^{k^r). This integral determines 
the geometric coupling of a triangle in Fourier space with a triangle on the CMB sky. 

Again in the limit of collapsed triangles, we can approximate the above coupling integral 

as 

Ci^i^i^{ki,k2,k3) ji^ikiro) j r'^drji^{k2r)ji^{k3r), (46) 
and using the definition of the 5-function, 

CiM^iki, k2, ks) ~ jii(fcirz))^^^^^-2-^. (47) 

Here the slowly varying spherical Bessel function is evaluated at the surface of last scattering 
since the transfer function Al_^{ki) in Eq. ()45|) is peaked there. 
Substituting this result into Eq. (jl^ we find 

~ fNL J fc2rffciA;2'rfA;2j;,(A;ir^)A^A;i)A^^(fc2)Af^(fc2)[P(A:i)P(fc2) + eye.]. (48) 
This can be evaluated as 

- fNLCl''C-\ (49) 

where Cj'^ is Cf, or Cf depending on the values of j and k. 
Now the (S/N)"^ for a given collapsed triangle is 

!,2 f2 (r<'^,ir~'jM2 



<i.i. {CI + Nl ) {CI + Nl ) {CI + Nf^ ) 
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(50) 



As long we can measure the appropriate temperature or polarization fluctuations on the 
scale k, {Ci^^f/{Cl + iV/J(Cf^ + iVfJ ~ 1 and Eq. ^ becomes independent of h- After 
integrating I2 up to Imax) while keeping li <^ I2, we find 

(|)' « fNLCiKrt^^, (51) 

where 

t\ = , ' (52) 

is the cross-correlation coefficient and again we assumed that we could resolve the primary 
CMB anisotropics on scale l\. 

If i = T, then on the very large scales we are considering = 1 and we recover 

(^)^ oc fji^Al'^^^, (53) 

the result we found for our toy model, Eq. (jSHl)- If i = E, then ~ 0.9 on large scales for 
models without reionization and ~ 0.5 for models with a significant reionization optical 
depth so our conclusion still holds. 

This explains the results we found in our Silk damping toy model, Eq. It is 

important to remember that we can only detect the primordial non-Gaussianity when we 
can resolve the primordial anisotropics. If the observed CMB power spectrum for the I 
modes from some collapse triangle is dominated by instrument noise or power induced by 
gravitational lensing, then cancellation of the Silk damped power spectra in Eq. (|Hn|) will 
not occur and we will observe the type of saturation seen in Fig. El 



IV. CONCLUSION 



The wealth of recent observational data has allowed the kinematics of the standard cosmo- 
logical model to be rigorously tested, now we must turn to the theory of its initial conditions. 
Slow-roll inflation has become the standard scenario used to explain the initial conditions 
of cosmology. While slow-roll inflation makes several predictions, the Gaussianity of the 
underlying curvature fluctuations may be the most robust and therefore should be tested. 
Since the CMB contains additional information in its polarization patterns, we calculated 
the increase in the Signal-to-Noise (S/N) ratio of the optimal cubic estimator when polariza- 
tion information is included. The improvement in WMAP is small, from /jvl ~ 13.3 for just 
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temperature to 10.9 for both T and E. Since WMAP is too noisy to measure E with high 
S/N we should not expect a large change. For Planck the improvement is from /atl ~ 4.7 
to 2.9 and we find that an Ideal experiment, with no instrument noise and infinitesimal 
beam width, improves from /jvl ~ 3.5 to 1.6. With an Ideal experiment cosmic variance 
limited up to / ~ 3000 it might be possible to observe the three-point function produced by 
non-linearites in General Relativity. 

We also explored how the four-point function induced by gravitational lensing would 
degrade our estimator's S/N. For WMAP there is very little effect, however the constraints 
from Planck can be reduced by 25%. For the next generation of CMB experiments which 
will measure / > 1000 with good sensitivity, the estimator used to constrain the primordial 
non-Gaussianity should be derived including the effects of gravitational lensing. 

The scaling with Imax was shown to be related to the total number of observed indepen- 
dent pixels on the sky, implying that (S/N)^ oc fskylmax^^i^max/lmin)- K the underlying 
distribution was Poisson the logarithim would be absent implying that the total {S/N^ 
simply scales as Npi^ — fsky^max- The functional dependence on Imax of our full calculation 
using the radiation transfer functions produced by CMBFAST agreed quite well with the 
prediction of our toy model. We showed that Silk damping did not reduce the signal avail- 
able from small scales appreciably. By using a toy model it was shown that this perhaps 
unexpected result is caused by the contribution from collapsed triangles. 

Our ability to constrain primordial non-Gaussianity on small scales crucially depends on 
our ability to measure the primordial anisotropics on those scales. While we focused on 
gravitational lensing in this paper, there are many other mechanisms (thermal Sunyaev- 
Zeldovich (SZ), kinetic SZ, Ostriker-Vishiniac, etc.. ) that produce additional Gaussian 
and non-Gaussian CMB fluctuations on these extremely small scales. The influence of these 
mechanism on our estimator is difficult to determine because of the non-linear physics in- 
volved and their highly non-Gaussian nature. Futher work is needed in order to understand 
how these effects will change the conclusions of this paper. These effects further strengthen 
the case for measuring the CMB polarization on small angular scales, as the above secon- 
daries are not expected to be significantly polarized. 
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